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We study the formation and detection at the next linear e^e~ collider of bound states of level-1 
quark Kaluza-Klein excitations Bkk within a scenario of universal extra-dimensions (UED). The 
interactions of such Kaluza-Klein excitations are modeled by an driven Coulomb potential. In 
order to obtain the threshold cross-section, we employ the Green function method which is known 
to properly describe the peaks below threshold and to yield a net increase in the continuum region 
(above threshold) relative to the naive Born cross-section. We study such effect at different values 
of the scale {R~^) of the extra-dimensions with an explicit calculation of the mass spectrum as given 
by radiative corrections. The overall effect is roughly 2.7 at = 300 GeV and goes down to 2.2 
at = 1000 GeV and a relatively large number of events is expected from N^vents ~ 2.5 x 10"* at 
R^^ — 300 GeV down to N events ~ 10^ at _R~^ — 1000 GeV at the anticipated annual integrated 
luminosity of Lo = 100 fb"'^. We finally discuss some potentially observable signatures such as the 
multilepton channels 2j + 2£ + f and 2j + 41 + for which we estimate statistical significance > 2 
for R-^ up to 600 ~ 700 GeV. 
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I. INTRODUCTION 

It is well known that as early as 1921 Theodore Kahiza 
proposed a theory that was intended to unify gravity and 
electromagnetism by considering a space-time with one 
extra space-like dimension [1]. A few years later Oscar 
Klein proposed that the extra space dimension (the fifth 
dimension) is in reahty compactified around a circle of 
very small radius [2]. These revolutionary ideas have 
thereafter been ignored for quite some time. However 
recent developments in the field of string theory have 
suggested again the possibility that the number of space 
time dimensions is actually different from D = A (indeed 
string theory models require D = II, i.e. seven addi- 
tional dimensions). In 1990 it was realized [3] that string 
theory motivates scenarios in which the size of the extra 
dimensions could be as large as i? « 10~^^ cm (corre- 
sponding roughly to electro-weak energy scale (~ TeV) 
contrary to naive expectations which relate them to a 
scale of the order of the Planck length Lp 10 ^'^^ cm 
(corresponding to the Planck mass Mp = ^ficjG w 10^^ 
GeV). See also [4]. 

Subsequently two approaches have been developed to 
discuss the observable effects of these, as yet, hypothet- 
ical extra dimensions. One possibility is to assume that 
the extra space-like dimensions are flat and compactified 
to a "small" radius. This is the so called ADD model [5] 
where only the gravitational interaction is assumed to 
propagate in the extra-dimension. A second possibility 
is contemplated in the Randall-Sundrum type of models 
where the extra dimensions do have curvature and are 
embedded in a warped geometry [6, 7]. 

Universal extra-dimensional models were introduced in 
ref. [8] and are characterized, as opposed to the ADD 
model, by the fact that all particles of the Standard 
Model (SM) are allowed to propagate in the (flat) ex- 



tra space dimensions, the so called hulk. Here to each SM 
particle X^^^ corresponds in this model a tower of Kaluza- 
Klein states X^") (KK-excitations), whose masses are re- 
lated to the size of the compact extra dimension intro- 
duced and the mass of the SM particle via the relation 
^ ~'~ 1^- All important aspect of the 

UED model is that it provides a viable candidate to the 
Cold Dark Matter. This would be the lightest KK par- 
ticle (LKP) which typically is the level 1 photon. Many 
aspects of the phenomenology of these KK excitations 
have been discussed in the literature. For reviews see 
ref. [9-12]. In particular KK production has been con- 
sidered both at the Cern large hadron collider (LHC) and 
at the next linear collider (ILC). Direct searches of KK 
level excitations at collider experiments give a current 
bound on the scale of the extra-dimension of the order 
R-^ > 300 GeV. See for example ref. [13]. At the Fermi- 
lab Tevatron it will be possible to test compactification 
scales up to ^ 500 GeV at least within some partic- 
ular scenario [14-16]. 

Lower bounds on the compactification radius arise also 
from analysis of electro-weak precision measurements 
performed at the Z pole (LEP II). An important fea- 
ture of these type of constraints is their dependence on 
the Higgs mass. A recent refined analysis [17] taking into 
account sub-leading contributions from the new physics 
as well as two-loop corrections to the standard model p 
parameter finds that R~^ > 600 GeV for a light Higgs 
mass (m/f = 115 GeV) and a top quark mass mt — 173 
GeV at 90% confidence level (C.L.). Only assuming a 
larger value of the Higgs mass the bound is considerably 
weakened down to i?"^ > 300 GeV for mn = 600 GeV, 
thus keeping the model within the reach of the Tevatron 
run II. The finding of this precision analysis are in qual- 
itative agreement with previous results [18], but are at 
variance with the conclusions of a recent paper [19] where 
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an analysis of LEP data including data from above the 
Z pole and two loop electro-weak corrections to the Ap 
parameter pointed to R-^ > 800 (at 95% C.L.). 

Equally important turn out to be the lower bounds 
from the inclusive radiative decay B — > X^j. It has 
been shown in ref. [20] that a refined analysis includ- 
ing in addition to the leading order contribution from 
the extra-dimensional KK states, the known ncxt-to- 
next-to-leading order correction in the Standard model 
(SM) gives a lower bound on the compactification radius 
> 600 GeV at 95% confidence level (CL) and inde- 
pendent of the Higgs mass. 

In this work we study the formation, production and 
possible detection of bound states of Kaluza-Klein n ~ 1 
excitations at e+e~ collisions. The production of bound 
states of KK excitation has been the object of some previ- 
ous work [21]. As compared to ref. [21] where the bound 
states production rates have been estimated by using a 
Breit-Wigner approximation, our study makes use of the 
method of the Green function in order to estimate the 
bound state contribution at the threshold cross-section, 
an effect which can be as large as a factor of three when 
considering strongly interacting particles. In describing 
the interactions that allow the formation of level- 1 KK 
bound states we assume that the level- 1 KK quark excita- 
tions interact via an driven Coulomb potential. This 
allows the use of analytic expressions for the Green func- 
tion of the Coulomb problem but it should be kept in 
mind that the results and conclusions about formation 
and decay of the bound state depend on this assump- 
tion. This method has also been recently used by the 
present authors in a study of sleptonium bound states 
within a slcpton co-next to lightest supersymmetric par- 
ticle (slepton co-NLSP) scenario of gauge mediated sym- 



metry breaking (GMSB) [22]. 

The plan of the paper is as follows. Section II briefly 
describes the UED model taken as a reference scenario. 
Section III discusses the formation criteria and shows 
that the bound states of KK level- 1 excitations do indeed 
form. Section IV describes the Green function method 
for the bound states providing an analytic formula for 
the Born production cross section. The threshold cross 
section for the bound state is studied for several values 
of the scale of the extra dimension R^^. Section V dis- 
cusses the possible decays of the bound states. Finally 
in Section VI we discuss the possible observation of the 
KK bound states at the e+e" linear collider pointing 
to three possible signatures whose standard model back- 
ground arc also considered providing an estimate of the 
statistical significance. In section VII we present the con- 
clusions. 



II. UNIVERSAL EXTRA DIMENSIONS 

The UED model is constructed considering the Stan- 
dard Model in a space time of 4 -I- D dimensions, and 
assuming that all SM particles are allowed to propagate 
in the extra dimensions which typically are assumed to 
be compactified to a radius R. In the following we fol- 
low strictly the notation of ref. [9] . We indicate the usual 
four dimensional coordinates as a;'', /i = 0, 1, 2, 3 and with 
a = 1, - ■ ■ D the extra space dimensions. The effective 
four-dimensional Lagrangian is then obtained by dimen- 
sional reduction, i.e. by integrating the A+D dimensional 
SM Lagrangian over the D extra space dimensions. Thus 
one has: 



+t (Q, u, d, L, e) {x^, y"") {T^D^ + t^+'^Da+a) (Q, u, d, L, e) (x^ , y'^) + 
Qix^y'') (Xu li(a;^y")^a2i^*(x^,^/'•)-f Ad y'^)iJ(a;^, y'^)) + H.c.l + 



(1) 



In the above Eq. 1 Ff- are the gauge field strength ten- 
sors of SM gauge group SU{2) x U{1) x SU{3) and 
are the gauge coupling constants in (4 -I- D)-dimensions 
which have dimension of (mass)^^^^ as well as the 
Yukawa couplings Xu,d.e- -D^ = d/dx^^ — and 
Da = d/dy'^ — A^+a are the covariant derivative and 



Aa = -iY,\=i9hTl{AAYi are gauge fields. Q,L are 
the SU (2) doublets, while it, d, e are the singlet (4 -I- D)- 
dimensional fcrmion fields. T^,A = 0, - • • (3 + £>) are 
(4+D)-dimensional gamma matrices satisfying the anti- 
commutation relations {r"^,r^} = 2g^^ . In the fol- 
lowing we shall deal with the simplest case of only one 
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extra-dimension {D = 1) ^. 

In order to extract the four-dimensional effective tlie- 
ory one needs to specify liow tfie extra-dimensions arc 
compactificd. It is found that in order to reproduce chiral 
fcrmions in 4-dimcnsion (the SM fermions) one is forced 
to assume an orbifold compactification structure which 
depends on the number of extra dimensions. For D-odd 
(e.g. D = 1) one chooses an S^/Z2 orbifold structure 
with Z2 being the reflection symmetry y — > —y. One as- 
sumes that the gauge fields Afj, and the Higgs boson H 
are even under the y — > —y transformation, while the ^5 
is assumed odd. This results in a Fourier series expan- 
sion of the fields which defines the zero modes (that corre- 
spond to the SM particles) and the level-n KK-excitation 
(coefficients of the expansions). 



if(a;^2/) 
A^{x^',y) 
A5 



00 



= U°(^,) + V2^^;:(x,)cos(^) 




(2) 



As already anticipated, difficulties arise when trying to 
construct chiral fermion fields in more than four dimen- 
sions. This is because in 5-dimensions, for example, it 
is not possible to construct the equivalent of the 75 ma- 
trix and bi-linear quantities like ipj^j^ip are not invariant 
under 5— dimensional Lorcntz transformations and there- 
fore they cannot appear in the 5D-Lagrangian. This ul- 
timately implies that for each standard model field one 
must introduce two 5D fermion fields whose zeroth or- 
der modes combine to give the AD chiral fermion. This 
however leaves some extra massless degrees of freedom at 
the zero level which can only be eliminated by formulat- 
ing the theory on an orbifold [10, 23]. The 5-dimensional 
fermion field is thus expanded as: 



vI/(x^y) 



1 



yj^^ixn + vi"^(x'')cos 



.V'^")(x^)sin(^ 



(3) 



These relations are however modified by radiative cor- 
rections which turn out to be cut-off dependent. These 
radiative corrections arise from loop diagrams traversing 
the extra-dimension [24] (bulk loops) and from kinetic 
terms localized on the brane which appear on the orbi- 
fold structure. 



i5(m^(„)) 

(S(mQ(„)) 
(5(m„(„) ) 
(5(mrf(„)) 
5(mi(„,) 
(5(me(„)) 
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(5) 



Here C,(z) is the Riemann zeta function, C(3) w 1.2020, 
and A is the cutoff scale of the theory. It correspond 
to the energy scale at which the effective 5-dimensional 
theory will break down, that is where the 5-dimensional 
couplings become strong and the theory is no longer per- 
turbative. A is the only additional parameter of the UED 
model beside the size of the extra dimension R. It can be 
estimated requiring that that loop expansion parameters 
remain perturbative. It has been found that the SU(3) 
interaction becomes non perturbative before the other 
gauge interactions for values of A > 10i?~^. So the par- 
ticle spectrum is typically computed with the above Eq. 5 
taking Ai? = 5, 10, 20. In the above expressions the brane 
kinetic terms are those dependent on the cutoff scale A. 
It should be noted that for KK scalars and spin-1 bosons 
the corrections in Eq. 5 simply add to Eq. 4, while the 
corrections for the fermion masses are introduced via the 
replacement ri? / R^ — >• {n/ R + Sm'^^^y . 



Performing the integration over the extra space dimen- 
sion the derivatives with respect to y will bring about 
mass terms that scale with the compactification radius 
R. Every level-n KK excitation A(„) acquires, in addi- 
tion to the SM mass (level 0), obtained via the Higgs 
mechanism, a new term: 

7^2 

TO^(„) = TO^(O) + -J, (4) 



^ When D=l, one ean choose = 7,1 and r4 = 175 



III. uiui BOUND STATE FORMATION 

In this section we shall review the possible creation of a 
bound state of the level- 1 Xi^T-excitation of the w-quark, 
i.e. a bound state uiUi. The interaction among two 
Kaluza-Klein excitations are driven by the QCD interac- 
tion, thus bearing no differences with respect to the Stan- 
dard Model; the strength of the interaction is given by 
as computed at a suitable scale [25-27]. We shall adopt 
the same formation criterion stated there, namely that 
the formation occurs only if the level splitting depending 
upon the relevant interaction existing among constituent 
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i?"' (GeV) 


KK mass (GeV) 




State mass M (GeV) 


Sis (GeV) 


A£(2P - IS) (GeV) 


300 


358.54 


0.136 


714.06 


2.937 


2.203 


400 


478.05 


0.131 


952.57 


3.627 


2.720 


500 


597.56 


0.127 


1190.92 


4.279 


3.209 


600 


717.08 


0.124 


1429.30 


4.903 


3.677 


700 


836.60 


0.122 


1667.69 


5.505 


4.128 


800 


956.11 


0.120 


1906.11 


6.089 


4.567 


900 


1075.62 


0.118 


2144.54 


6.658 


4.993 


1000 


1195.14 


0.116 


2382.98 


7.214 


5.411 



TABLE I: Results of Coulombic model for the bound state of the level-1 iso-doublet f/i quark. The strong coupling Qs is 
computed at the scale Q = r^^, where tb = 3/(2mas) is the Bohr's radius. For each mass value m the scale Q = r~j^^ 
depending itself on as must be solved numerically from the equation Q = (2/3) mas{Q). 



particles is larger than the natural width of the would-be 
bound state. This translates into the formation require- 
ment 



AE2P-1S > r 



(6) 



where Ai?2P-is ~ E2P — Eis and F is the width of the 
would-be bound state. The latter is twice the width of 
the single KK quark, F = 2Tkk, as each KK quark 
could decay in a manner independent from the other. 

We shall stress that F bears no resembling to the to- 
tal decay width of the resonance, as it only includes the 
single KK decay mode and not the annihilation modes 
discussed later in sec. V. It represents the minimal en- 
ergy level spread needed for bound state formation, which 
allows for the separation among the fundamental and 
the first excited state. If, and only if, the bound state 
is formed, then it is possible to discuss its annihilation 
widths as described in sec. V. In our model V{x) is given 
by a Coulombic potential 



V{r) 



4as 
3r 



(7) 



with r = I a; I, and where as is the usual QCD coupling 
constant which has been taken at a suitable scale as de- 
scribed in [25, 26]. This model has proved to be reliable 
because of high mass values involved in the problem, and 
gives the great advantage of having full analytical results. 
We are thus able to compute its energy levels given by 
the expression 



4 raaj. 



(8) 



and the separation of the first two energy levels is given 
by 



2P-1S 



(9) 



The scale at which is evaluated is given by the inverse 
of Bohr's radius rp = 'i/{2mas), the average distance of 
the constituents of the bound state. Therefore it is found 



by solving numerically the equation Q = (2/3) m a,, (Q) 
and it is of order 0(10 GeV) for m w 300 GcV and 
0(100 GeV) fr m w 1200 (GeV). The corresponding val- 
ues of as are given in Table I. The mass of the nth bound 
state is given by the expression: 

M„ = 2m + e„ (10) 

where m is the mass of the constituent ui quark and £"„ 
is given by (8). The wavefunction at the origin, which 
will be needed in order to compute decay widths, for this 
particular model is given by the expression 



1^(0)1^ 



1 



(11) 



The obtained results are given in Table I. 

We observe that the bound state energies are of the 
order of the GeV for this range of KK mass, and that 
the spreading of the first two bound states raise linearly 
with m. 

In order to determine whether the bound state will 
be formed we shall apply the criterion given in eq.( 6). 
The KK-q\xaA decay widths have been already com- 
puted in [21], where it has been shown that their values 
are at most of the order of 100 MeV, one order of mag- 
nitude less than the energy splittings. In this scenario 
the eq. (6) requirement is always fulfilled, and the bound 
state is formed for KK-q\xaik masses in this investigation 
range. 

IV. GREEN FUNCTION 

In order to describe the cross-section of a KK bound 
state in the threshold region we shall use the method of 
the Green function. We bricfiy review here the essential 
features of the mechanism, and refer the reader to the 
literature for further details we [28]. We start from the 
Schrodinger equation which describes the bound state by 
means of a suitable potential V{x), 



(12) 
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where £ is the energy eigenvalue of the bound state. The 
threshold cross-section of this bound state is then pro- 
portional to the imaginary part of the S wave Green func- 
tion of this Schrodinger equation, Qis{x^ y, E), where the 
two constituent particles are sited in x,y and E is the 
energy offset from the threshold (not to be confused with 



iH-E)gisix,y,E)=5{x-y) 



(13) 



By means of the substitution E ^ E + iT wc take into 
account the finite width of the state. 

The cross-section is thus proportional to the expres- 
sion 



Im 



(14) 



x=0,y=0 



the derivative of eq. (14) has a simple expression, as we 
have 



9^is(0,0,£;) . (15) 



a;=0,y=0 



The complete expression for the IS" Green function of 
our problem function of energy from threshold is 
given with a slight change of notation by [30]: 



GisiQ^G^E + iT) 



m 
An 



-2A 



k 



log 



.2A 

V'(i - i^) + 27 - 1 



(16) 



where k = ^ —m(E + iP), A = 2Q;sm/3 and the wave 
number \s v ~ Xjk; E = ^/s — 2m. The ip is the 
logarithmic derivative of Euler's Gamma function T{x), 
7 ~ 0.57721 is Euler's constant and is an auxiliary pa- 
rameter coming out from a dimensional regularization, 
the factorization scale, that cancels out in the determi- 
nation of physical observables. 

The final expression for the production cross-section 
of a KK bound state is thus given by 

cr(7Ti, E, r, Us) = ^^(TB Im [Qis] (17) 

where ag is the Born expression of the cross-section [31] 
The process e+e~ UiUi proceeds through the anni- 
hilation into the standard model (level-0) gauge bosons 
7 and Z but in principle one should also consider the 
contribution of the level-2 gauge bosons 7(2) and ^(2)- 
Especially so in our case of threshold production of the 
pair uiUi. Indeed in this case m w 1/-R, and ^/s = 
2m + E K. 2 /R + E and since m^^ fa 2/R when pro- 
ducing at threshold the uiUi pair we would be close 
to the 72 and Z2 resonances. However as discussed in 
section II the mass spectrum is modified by the radia- 
tive corrections. We have verified that over the region 
of parameter space 300 GeV < R'^ < 1000 GcV and 
2 < AR < 70 the pair production threshold 2m„j is 



always larger than rUj^ , and thus these resonances 
should in principle be included in the calculation. We 
have also verified, cross checking our calculation with the 
output of a CalcHEP [29, 32] session, that the numerical 
impact of these diagrams is completely negligible. Their 
contribution turns out to be five orders of magnitude 
smaller than that of the SM gauge bosons 'y, Z. The an- 
alytic formula of the Born pair production cross section 
e+e~ 7*,Z* UiUi can be deduced for example 
from those of heavy quark (tt) [33] taking into account 
the fact the the level-1 KK quarks are vector-like i.e. 
their coupling to the Z is of the 7^ type and has no axial 
component. Following the notation of [33] the amplitude 
is written as: 

M ^ v{k2)i^{av + bvf)u{ki) X 

Dv[s, Mvj 

where Dy{s) = s — My + iMz^z is the gauge boson 
propagator factor, ay, by arc the [standard model) cou- 
pling coefficients of the electron to the gauge bosons while 
Ay, By are the coupling coefficients of the level-1 KK Ui 
quark to the gauge bosons. These electron coefficients 
are: = —1, = 0^ az = —1/4-1- sin^ 9w, bz = +1/4, 
with 9w the Weinberg angle of the SU{2) (g) U{1) gauge 
theory. The Ui coefficients are: = +2/3, Az = 
1/2 - {2/3)sm^ew [9], = Bz ^ (recall that Ui is 
vector-like). Finally = e and gz ~ e / (cos 9w sin 9w) 
with e the electronic charge. The final expression of the 
Born pair production cross section is: 



c^B 



c^z + o-^z 



az = 47ra" 



/3(s + 2m2) 



r 2 212 



-Sttq; 



(a| + b\)A\ / 3(g + 2m^) 

\Dz{s)\^ 
^azA^Az I3 {s + 2m^) 
\Dz{sW 



(1 



M; 



(18) 
(19) 

(20) 
-) (21) 



where cw = cos Ow and sy/ = sin6'H' and a — / (47r) is 
the QED fine structure constant. From Eq. (16) one can 
readily see the behavior of the cross-section (17) for large 
E is given by k. The finite width of the state has been 
taken into account by the substitution E ^ E -\- iV ^ and 
this position makes a great quantitative difference below 
threshold. When computed for positive energy offset the 
variation of F makes essentially no difference for the re- 
sulting cross-section. 

In this work we shall concentrate on the continuum 
region of the cross-section, namely E > 0. The region 
below threshold, E < 0, has been already discussed in 
detail in ref. [21], where the authors presented an analysis 
of both the positions and the widths of the peaks using 
a Breit-Wigner description. In this respect the Green 
function approach does not carry substantial differences 
relative to the Breit-Wigner one. Indeed the position of 
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a) m[/, =358.541 GeV, R-^=300G&\ b) mj;, -478.054 GeV, R-^=400Ge\ 




02468 10 02468 10 

E (GeV) E (GeV) 



FIG. 1: Production cross sections of level-1 KK doublet quarlc bound states UiUi as a function of the energy offset from 
threshold (yG = 2mui + E), for values of the scale of the extra-dimension in the range [300 -i- 600] GeV and a total width 
of r = 0.5 GeV. The continuous line is the Green Function result, the dotted one is the Born approximation given by our 
analytical formula (Eq. 18). The full circles represent the Born cross section from the CalcHEP [29] numerical session inlcuding 
also the annihilation diagrams of 72 and Z2 whose contribution is however completely negligible. The numerical results from 
CalcHEP are in complete agreement with our analytical formula in Eq. 18. The cut-off scale A, at which perturbative expansions 
break down, has been fixed so that AR — 20. 



the poles and the broadening of the peaks are the same 
due to the presence, inside the ip function of Eq. (16), 
of terms which include the binding energy En of Eq. (8) 
and the decay width T. For values of E close to En, 
the argument of the ip function inside Eq. (16), namely 
(l — v), approaches a negative integer, simple pole of the 
function in the complex plane, while the presence of the 
r determines the width of the peak centered in En- 

In Fig. 1 we show the cross-section for a range of 
the value of the scale of the extra dimension, R^^ = 
300 — 600 GeV, while Fig. 2 provides the same plots are 
shown the range R'^ = 700 - 1000 GeV. In both figures 
the value of the other parameter is fixed at AR = 20. 
This parameter enters our calculations only when com- 
puting the mass spectrum through the logarithmic terms 



in Eq. 5. We thus provide a quantitative study of the 
effect of the formation of bound states of the level-1 KK 
quarks with respect to the parameter of the model (R^^). 
The results are less sensitive to the other parameter {AR) 
which only enters through the logarithmic factors in the 
radiative correction terms in the mass spectrum of the 
model. In Figs. 1& 2 we have fixed AR = 20 and varied 
R~^ computing the corresponding values of the level-1 
KK quark mass, and assuming the energy of the collider 
being fixed at y/s = 2mu^ +E, E being the energy offset 
from the threshold. We have used a value of F = 0.5 GeV 
for illustrative purpose, compatible with the formation of 
bound state. Different choices of F by even two orders of 
magnitude smaller will not make a visible difference on 
the figures. 
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a) mj/ =836.6 GeV, 7?-i=700GeV 



b) mj/ =956.11 GeV, R-^=SOOGqY 




2 4 6 8 10 

E (GeV) 

c) mt/j =1075.62 GeV, i?-i=900GeV 





4 6 8 10 

E (GeV) 

d) mt; =1195.14 GeV, i?-i=1000GeV 




FIG. 2: Production cross sections of level-1 KK doublet quark bound states UiUi as a function of the energy offset from 
threshold (y^ — 2mui + E), for values of the scale of the extra-dimension in the range [700 -ir 1000] GeV and a total 
width of r = 0.5 GeV. The continuous line is the Green Function result, the dotted one is the Born approximation given by our 
analytical formula (Eq. 18). The full circles represent the Born Cross section from the CalcHEP [29] numerical session inlcuding 
also the annihilation diagrams of 72 and Z2 whose contribution is however completely negligible. The numerical results from 
CalcHEP are in complete agreement with our analytical formula in Eq. 18. As in Fig. 1 we have kept fixed AR — 20. 



One can observe that the cross-section obtained by 
the Green function has a behavior like Ve for small 
energy offset. The cross-section decreases with increas- 
ing mass: for E ^ 10 GcV its value is about 250 fb 
at R^^ = 300 GcV; goes down to approximately 70 fb 
at R^^ = 500 GcV. Finally it approaches 13 fb at 
R~^ = 1000 GeV as can be seen form Fig. 2. The 
Green function cross section is larger than the Born cross- 
section by a factor that ranges from 2.7 (at R^^ = 300 
GeV) down to 2.2 (at R'^ = 1000 GeV) at the same en- 
ergy offset value {E = 10 GeV). This result is due to the 
fact that the Green function method takes into account 
the existing interaction among constituent particles, and 
the contribution of binding energies accumulate towards 
the e = level, thus substantially contributing to the 
continuum region as well. 



An important consideration is in order here. The 
results for the bound state cross-section given depend 
solely on the coupling constant as{rB) and the mass of 
the KK excitation, thus they are universal to bound 
states made out of other flavors of KK quarks. This 
however docs not apply to other kinds of KK excita- 
tions bound states, like for instance bound states of KK- 
leptons. In this case we have a = uqed for coupling con- 
stant, much weaker at this scale than the strong coupling 
constant aa{rB)- The QED coupling constant would 
lead, not only to lower values for the absolute production 
cross-section, but will also reduce drastically the main ef- 
fect being discussed here, i.e. the enhancement, above 
threshold, due to the bound state interaction relative 
to the Born cross-section. The threshold cross-section 
would be only a few percent larger than the Born cross- 
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FIG. 3: Solid line: decay width (keV) of the pseudoscalar 
uiUi bound state to two photons as a function of KK exci- 
tation mass. Dashed line: decay width (MeV) of the pseu- 
doscalar ui'Ui bound state to two gluons as a function of the 
KK mass. 



section. It would not clearly be an effect as large as that 
shown in figs (1), (2) which turns out to be quite strik- 
ing, i.e. the threshold bound-state cross section is about 
three times as large as the Born result. 



V. uiui DECAY WIDTHS 

The KK bound states we discuss here arc the pseu- 
doscalar -^Sq and the vector one ■^S'l . For the pseudoscalar 
state the decay channels arc into two photons or two glu- 
ons for which the following Born level expressions hold 
(see for instance [26]): 



^bCSo 77) = qfa 



4 248^I^(0)P 



M2 



and 



^bCSq gg) 



(22) 



(23) 



Here qi is the charge of the constituent quark of the 
bound state, while M and IV'IO)!^ are given by (10) 
and (11) respectively. 

The QCD radiative correction [34], which is the same 
in the two cases, lead to the following one-loop width: 



r? 



TT 



20 



The results obtained for the two decays ^Sq 
^Sq ^ gg are shown in Fig. 3. 



(24) 



77 and 



> 400 




400 600 800 1000 
Mkk (GeV) 



1200 



FIG. 4: Decay widths of the uiUi vector bound state to two 
charged particles and into three gluons as a function of KK 
mass. Here we have considered all possible e.m. decay chan- 
nels. 



For the vector case '^S'l the relevant decay channels are 
the one in charged pairs and the one into three gluons, 
for which one has 



2 2 2l6^|V(0)|^ 
q^q^a 



M2 



and 



TBCSi^ggg) 



3l607r|^.(0)|2 



81Af2 



(25) 



(26) 



The charge of the final state charged particle is given 
by qf. The QCD radiative corrections [34] modify these 
expressions into 



and 



q+qj)^rBCSi^q+qj)[l~ 



16 



(27) 



ggg)^TBeSi 

TT 



X 



1 



14 



.9.9.9) X 
27 



(1.161 + log(2)) 



(28) 



Observe that the as that appears in the perturbative 
corrections has to be computed at a scale of the order 
of 2m. It is thus different from the as occurring in the 
expression of the wavefunction at the origin given by (11), 
the latter being computed at the scale of the inverse of 
Bohr radius. 

The two decays of the vector state are shown together 
in Fig. 4. 
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We observe that only the pseudoscalar hadronic decay 
is in the MeV range and raises approximately linearly 
with KK mass. The ^Sq photonic decay and ^Si decays 
are smaller by almost two orders of magnitude for the 
considered KK mass range. For the pseudoscalar case 
the hadronic is the dominant decay by far, while in the 
vector case the decay into charged particles, when tak- 
ing into account all possible processes as seen in Fig. 4 
overtakes the hadronic decays. 

Other electro-weak decay channels are negligible. 
Those arc proportional to , thus their ratio to gluonic 
decays is suppressed by (a/asY, at least by two orders 
of magnitude. 

For most scenarios depending upon the values of A and 
R [21] single quark decay becomes the dominant decay 
channel for the bound state. 

Any two-body decay width of the bound state is pro- 
portional to 3^1-0(0)1^/^/^, where g is the relevant cou- 
pling to the decay particles. Thus any electro- weak width 
is in the keV range, as previously seen. This result is true 
in general for any two-body decay process, the only no- 
table exception being the hadronic decay of (23), as the 
g coupling this time is rather large, being equal to . In 
this case the value is in the MeV range, as seen in Fig. 3. 

Three-body decays are further suppressed with respect 
to previous formula by another power in g and phase- 
space reduction, resorting again in the keV range of en- 
ergies. 

From [21] one sees that in most cases single quark de- 
cays (SQD) are by far the most important decay chan- 
nels of the bound state, to the order of hundreds of MeV, 
while as discussed above bound state decays are essen- 
tially negligible. Moreover a comparison of those SQD 
widths with the results of Table I through eq. (6) shows 
that for the considered mass range of KK there is for- 
mation of the bound state. 



VI. DETECTION 

As we have previously seen, for large values 
{R~^ > 300 GcV) SQD is the dominant decay channel 
for a KK bound state, thus leading a to a dominant 
signature consisting of two monochromatic quarks plus 
missing energy. Other interesting signatures that could 
be considered, for example the three jet production due 
to the bound state decay into three gluons ^Si ^ ggg 
discussed in the previous section, are clearly subdomi- 
nant given the fact that B{^Si ggg) sa 10~^ ^. Such 
signatures, in addition would have to be confronted with 
important QCD backgrounds. Therefore in the follow- 



^ The value B{^Si — ^ ggg) ~ 10~^ is easily obtained assuming a 
SQD width of the order of a few hundreds MeV and combining 
this with the results of the partial widths given for example in 
Fig. 4 (or by using Eq. 25 and Eq. 26). 



ing discussion we concentrate on the dominant channels 
given by the single quark decay whose branching frac- 
tions, on the contrary, can be as high as 65% and 98% 
and involve missing energy in the final state. Follow- 
ing [28] we limit our analysis to the region above thresh- 
old, i.e. E > 0. The region below threshold, E < 0, 
is characterized by peaks in the cross section for values 
of E equal to binding energies of the bound states. The 
width of those peaks are given by the decay width of the 
bound state, which are at most of the order of the MeV 
for the SQD and much less, of the order of the keV, for 
other annihilation decay modes, as discussed in sect. V. 

From Eq. (8) we can estimate the separation of the 
various peaks below threshold, which tend to merge when 
they accumulate, that is for n such that 



—ma, 
9 



1 



1 



(n + l)2 



(29) 



in this manner we estimate that the last resolved peak 
has a quantum number n that satisfies 



2n + l 9r _ r 

n2(n + l)2 4ma2 Ei 



< 1 



(30) 



from the values Table I and using a width value of the 
order of 20 MeV there are only around 5 peaks left before 
merging. 

Because of ISR and beam energy spread, of the order of 
the GeV for a future linear collider, it is unclear whether 
it could be possible to resolve those peaks of keV magni- 
tude with this machine. The only potentially detectable 
peaks should be the ones belonging to a SQD, provided 
one has a scenario with widths of the order of the MeV. 

The situation above threshold changes drastically with 
respect to the "naive" Breit-Wigner estimate, as is 
clearly shown in Fig. 1 and Fig. 2 . A few GeV above 
threshold make for a factor of 3 of increase compared 
to the Born cross-section, allowing a clear distinction 
between the two cases. Assuming an annual integrated 
luminosity of Lq = 100 fb^^ and a scale of the extra- 
dimension R^^ ~ 300 GeV one finds around 2.5 x 10''^ 
events per year of two quark decay for a center of mass 
energy of 10 GeV above threshold (we adopt here the 
scenario for which the branching ratio of SQD is essen- 
tially 1). The number of events per year loses an order 
of mag nitude at R'^ = 700 GeV, that is about 3 x 10^, 
as could be inferred from Fig. 2. 

As we have already said the decay width of the KK 
bound state will be given by twice the decay of the single 
quark, as the SQD dominates, being of the order of up 
to hundreds of McV. For our uiui bound state there are 
two possible scenarios of decay pattern [24] . The first one 
concerns the iso-singlet ui^ for which the decay channel 
into Wi is forbidden while that into Zi is heavily sup- 
pressed B{uij^ — > ZiUor) sin^ 6'i « 10^^ 10^'^ and 
the dominant channel is given by — > uqrTIj with 
B{uij^ — >■ uojj7i) ~ 0.98 whose signature is a monochro- 
matic quark and missing energy of the KK photon, the 
latter being the LKP [24]. 
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For the iso-doublet uij^ the situation is more interest- 
ing, as more channels are available [24], notably uij^ — > 
doj^Wi with B{ui^ (Iq^Wi) ~ 0.65 and ui^ uq^Zi, 
with B{ui^ — > uq^Zi) « 0.33 while the branching ratio 
into 7i is negligible B{ui^ — >■ uo£,7i) ~ 0.02. The decay 
chain into Wi can follow the scheme: ui^ — > do^Wi 
c^Ol^oi^Il ~^ doj^iovoli with branching ratio given by: 

B{ui^ dojoi'oli) ~ B{ui^ do^Wi)B{Wi Iqvi) 
w 0.65 X 1/6 X 1 « 10"^ (31) 

and alternatively, the same final state could be reached by 
the scheme: ui^ -> do^Wi do^iiVoL ^Ot^o^'OiTi- 
As compared to the iso-singlet case, the result is a 
monochromatic quark, a lepton and missing energy in 
both cases. 

The decay into the Zi channel is ui^ — > uq^Zi — >■ 
uql^o^i ^ uo^^o^o7ij resulting in a monochromatic 
quark, two leptons and missing energy. The branching 
ratio of the above chain is: 

B{ui^ uot447i) ~ B{ui^ uo^Zi)B{Zi LqLi) 

xS(Li^47i) 

~ 7:X-xlw5x 10^2 (32 
3 6 

These leptonic decays of ui have much cleaner signatures 
than the hadronic ones allowing, in principle, for a better 
detection of the signal. 

In all cases we emphasize that the observable signal of 
the bound state production at the linear collider would 
be similar to that of the Born pair production except 
for the absolute value of the cross-section. In particu- 
lar, assuming for definiteness a linear collider operating 
around the threshold ^/s = 2m + E GeV we would have 
for R-^ = 300 GeV and E = IQ GeV, in the case of an 
iso-singlet bound state (or Born pair production of uiji): 

e+e" ^ 2jets + |; (33) 

with cross section: 

(T(e+e" 2jets + gb^^ x [B{uii{ -> wo7i)]^ 

« 173 fb. (34) 

We note that the (Jbk.ic for the iso-singlet ui has to be 
computed ex-novo and cannot be read from the values of 
Fig. 1 since it refers to the iso-doublet Ui. The singlet 
and doublet have, when including radiative corrections, 
different masses and the corresponding pair production 
threshold is therefore different. For the values of the scale 
of the extra-dimension R^^ = 300 GeV the mass of the ui 
iso-singlet is — 351.75 GeV (slightly lighter than the 
iso-doublet) and the corresponding Green Function cross- 
section at an energy offset of i? = 10 GeV is cre^^ = 181 
fb. At an e+e" collider this signal has a standard model 
background from ZZ production with one Z decaying to 
neutrinos and the other decaying hadronically. The cross 
section for ZZ boson production is « 244 fb at an energy 



offset of = 10 GeV from the relative thresholds. This 
provides the following estimate for the SM background 
at = 300 GeV: 

crsM(2jets + |:) « 244fb X 0.7 x 0.2 « 34fb (35) 

In the case of an iso-doublet bound state (or Born pair 
production of Ui) the Wi decay chain gives the signal: 

e+e" ^ 2jets + 2^ + |; (36) 

with cross section (see Fig. 1): 

aie+e- ^ 2j + 2£ + |^ = ctb^^ [S(uil ^ doJoi^oli)f 
w 275fb X (10~^)2 
= 2.75 fb (37) 

while the Zi decay chain gives rise to the signature: 

e+e" ^ 2jets + 4£ + |^ (38) 

with cross sections: 

aie+e- ^ 2} + M + = ae^^ [B{uil ^ uo^447i)]' 
= 275fb X (5 X IQ-^f 
« 0.69 fb (39) 

Triple gauge boson production, WWZ, ZZZ at a high 
energy linear collider has been studied in refs. [35, 36]. 
It has been found that these processes receive a sub- 
stantial enanchement in the higgs mass range 200 GeV 
< rriH < 600 GeV particularly the ZZZ channel. As 
these processes provide a source of standard model back- 
ground for our signal we estimate them both at a value 
of m/i = 120 GeV and at a value of ruh = 200 GeV 
for which the cross sections are enhanced. Production 
of WWZ can for instance give rise to the signature of 
2jets + 2£ + ^ via leptonic decay of the W gauge bosons 
and hadronic decay of the Z boson, while the ZZZ pro- 
duction can produce 2jets + A£ + f! via hadronic decay 
of one Z while the others decay leptonically with one 
of them to a pair of r which subsequently decay to ii^D 
{£ = e,iJ,). Estimates of the resulting cross sections are 
found using the CalcHEP [32] and CompHEP [37] soft- 
ware. We have verified agreement with previous results 
given in ref. [36] and for a Higgs mass of nih ~ 200 GcV 
we find, for R-^ = 300 GcV at = 2mu, + 10 w 724 
GeV: 

a(WWZ) « 72fb (40) 
a{ZZZ) « 7fb (41) 

We thus find at = 724 GeV within the standard 
model: 

asM{2i + 2e + ^ w 70 fb X (0.1)2 X 0.7 « 0.5 fb 

(Tsubi + M + f:) w 7 fb X (0.3)2 X 0.7 X (0.17)2 

« 1.2 X IQ-^ib (42) 

The 2jets-|-2£-|-^channel could be potentially contami- 
nated also from tt pair production cross section which at 
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R-^ (GeV) 


(GeV) 
(iso-singiet) 


2jcts + |; 


2jcts + 2^ + |; 
ruh = 120 GeV (200 GeV) 


2jets + 4^ + |; 
ruh = 120 GeV (200 GeV) 


300 


351.7 


121.8 


13.1 (12.9) 


8.3 (8.2) 


400 


469.0 


81.1 


8.2 (8.1) 


5.6 (5.6) 


500 


586.2 


58.7 


5.7 (5.6) 


4.2 (4.2) 


600 


703.5 


44.4 


4.1 (4.1) 


3.3 (3.3) 


700 


820.7 


35.3 


3.1 (3.1) 


2.7 (2.7) 


800 


938.0 


29.0 


2.4 (2.4) 


2.3 (2.3) 


900 


1055.2 


24.1 


1.9 (1.9) 


1.9 (1.9) 


1000 


1172.5 


20.8 


1.6 (1.6) 


1.7 (1.7) 



TABLE II: Estimate of the statistical significance 55* as defined in Eq. 43 and corresponding to the annual integrated luminosity 
I/O = 100 fb~^ for the three channels discussed in the text as a function of R^^ and ^/s = 2mui +E, assuming an energy offset 
of i5 = 10 GeV from the threshold. The physical threshold of the 2 jets + channel is different form that of the other two 
channels as it refers to the ui iso-singlet level-1 KK quark whose masses at various values of R~^ are given in column two and 
can be compared with the corresponding masses of the Ui state form Fig. 1 and Fig. 2. The values of the statistical significance 
for the two multilepton channels have been computed for two values of the Higgs mass nih = 120 (200) GeV. Again A_R = 20. 



such high energies is 0(300) fb [38]. Assuming the top 
quarks to decay with probability one to Wb and then 
the W gauge boson decay via the leptonic mode (with 
B{W — )■ ii/g) ss 0.1) would mimic the signal with a cross 
section (T5M(2jets + 2£ + fJ) ^ 3fb. However in this case 
we expect 6-tagging of the hadronic jets. Assuming an 
efficiency in 6-tagging of 60% we would get a contribution 
of 1.2 fb to the 2jets + 2£ + cross-section which has to 
be added to that in Eq. 42. This has been done in the 
calculation of the statistical significance of table II. 

We conclude providing an estimate of the statistical 
significance: 

SS = (43) 

of the three signals discussed above as related to an in- 
tegrated luminosity of Lq = 100 fb^^ {Ng is the number 
of signal events and is the number of background 
events). These estimates are given in table II and Fig. 5. 
Albeit quite encouraging (especially so the SS of the 
2jets -I- we should bear in mind that the actual obser- 
vation of these signals might be not be so easy from the 
experimental point of view. Indeed it is quite likely that 
in a framework of a quasi degenerate KK mass spectrum 
the jets will be typically quite soft and therefore difficult 
to detect. It is therefore customary to concentrate on 
the much cleaner multilepton signatures [24, 39]. Indeed 
a similar analysis to the one given here, but with a per- 
spective on signals arising at the Compact Linear Collider 
(CLIC), regarding the (Born) pair-production of level-1 
KK-leptons and level-1 KK- quarks is given in rcf. [39]. 

It is also well known that jets, multilepton and miss- 
ing energy signals are as well typical of supersymmetric 
models. Indeed detailed studies have already appeared in 
the literature regarding the possibility of distinguishing 
supersymmetric and universal extra dimension models at 
both the large hadron collider and linear collider: sec for 
example rcf. [40-42]. 
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300 400 500 600 700 800 900 1000 
^-'(GeV) 

FIG. 5: Statistical significance (SS) for the various channels 
discussed in the text as a function of the scale R~^ of the 
extra-dimensions. Note that the statistical significance of the 
2j + is scaled down by a factor of 10. 



However, in all cases, angular, invariant mass and/or 
missing energy distributions of the discussed signals 
would be identical to those obtained in the Born pair 
production. In our opinion further detailed analysis of 
the signals and of the possible SM backgrounds {and/or 
competing SUSY signals) goes beyond the scope of this 
study, whose main objective is to emphasize the dramatic 
increase of the bound state cross-section relative to the 
Born pair production. 



VII. CONCLUSIONS 

Within a universal extra-dimensional model we have 
considered the formation and decay of a bound state of 
level-1 quark Kaluza-Klein excitation and its consequent 



12 



detection at a linear e+e~ collider. Since itikk should 
be larger than at least 300 GeV we have used a model 
with a Coulombic potential. Admittedly this is a model 
assumption and it should be kept in mind that our re- 
sults are strictly valid only within this premise, which 
however has the advantage of providing full analytical 
expressions for the effect. Being a bound state we have 
used the Green function technique for the evaluation of 
its formation cross-section in the threshold region, which 
is more appropriate than the standard Breit-Wigner pic- 
ture as it takes into account the binding energy and the 
peaks of the higher level excitations that coalesce towards 
the threshold point. The net effect is a dramatic increase 
of the cross-section in the continuum region right of the 
threshold. This multiplicative factor is roughly 2.7 for 
300 GeV and drops down to 2.2 at i?"^ = 1000 
GcV. The Green function cross-section would allow more 
than « 10^ events per year even at R^^ = 400 GeV 
{mui ~ 478 GeV) for a suitable integrated luminosity of 
the e+e" linear collider {Lq = 100 fb^^). The number of 
events R'^ = 1000 GeV {mu, w 1200 GeV) would stih 
be « 10^ at the same integrated luminosity. 

The large difference among the two descriptions of the 



cross-section should also possibly help in the determina- 
tion of the correct model for such a heavy bound state 
outside the SM. 

Our analysis of the backgrounds to the final states sig- 
nals, though very simplified, indicates that the multi- 
lepton channels have a good statistical significance {SS > 
2) at least up to R'^ = 600 700 GeV, which cer- 
tainly warrants further detailed and dedicated studies of 
these channels and their backgrounds. The potentially 
large (one order of magnitude) estimated statistical sig- 
nificance of the 2j + fJ channel must be taken however 
with great caution because this type of signal may prove 
difficult to observe as it will be characterized by soft jets 
within the relatively degenerate mass spectrum of the 
extra-dimensional model. Further detailed studies are 
also needed for this channel. 
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